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Abstract 

Fixed-order QCD radiative corrections to the vector-boson and Higgs associated production 
channels, pp — > VH (V = W^,Z), at hadron colliders are well understood. We combine higher 
order perturbative QCD calculations with soft-gluon resummation of both threshold logarithms 
and logarithms which are important at low transverse momentum of the VH pair. We study the 
effects of both types of logarithms on the scale dependence of the total cross section and on vari- 
ous kinematic distributions. The next-to-next-to-next-to-leading logarithmic (NNNLL) resummed 
total cross sections at the LHC are almost identical to the fixed-order perturbative next-to-next-to- 
leading order (NNLO) rates, indicating the excellent convergence of the perturbative QCD series. 
Resummation of the VH transverse momentum (pt) spectrum provides reliable results for small 
values of px and suggests that implementing a jet-veto will significantly decrease the cross sections. 



I. INTRODUCTION 



The recent discovery of a Higgs-like particle has brought our understanding of 

electroweak symmetry breaking to a deeper level. Now it is imperative to study the detailed 
properties of this particle in the hope of finding any hints for new physics beyond the 
Standard Model (SM). An important Higgs production mechanism at hadron colliders is the 
associated production of a Higgs boson and a vector boson, VH (V = W ± , Z) [3|. At the 
Tevatron, the process qq' — > VH with the decay of the vector boson to leptons and of the 
Higgs to the bb and W + W~ channels has provided important sensitivity to a light Higgs 
boson ~\ 5|. At the LHC, the production rate for associated VH production is small, but 
with ~ 30 fb _1 a light Higgs in association with a W or Z can potentially be observed in 
the boosted regime via H — > bb p]. Reliable predictions are essential for the observation 
and study of the VVH couplings in this channel [7, 8]. 

The rate for associated VH production is perturbatively known to next-to-next-to-leading 
order (NNLO), i.e. 0(a 2 s ) 2I, Q. At next-to-leading order (NLO), the QCD corrections are 
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identical to those of the Drell-Yan process for an off-shell gauge boson, qq' — > V* 
At NNLO, however, the ZH process receives a small additional contribution from the qq 
initial state, gg — > ZH [9]. The NLO rates are available in the general purpose MCFM 14 1 
program, while the total rate can be found to NNLO using the VH@NNLO code [9J. 

Infrared finite results in higher-order QCD processes occur due to a cancellation of virtual 
and real soft divergences. The fixed-order calculation is reliable providing all of the scales 
are of the same order of magnitude. When the invariant mass Mvh of the final state 
particles WH or ZH approaches the center-of-mass energy of the colliding partons, there is 
less phase space available for real emission. While the infrared divergences will still cancel, 
large Sudakov logarithms will remain. These logarithms can spoil the convergence of the 
perturbative series and need to be resummed to all orders for reliable results in this threshold 
region 15|. Threshold corrections involve terms of the form a" log ^_^~ z ^ , which are large 



when z = My H /s ~ 1, where s is the partonic center-of-mass (cm.) energy-squared 
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20j. Similarly, large logarithms of the form a" log 



2n-l / M- 



~Pt,vh 



I can also occur when the 



VH system is produced with small transverse momentum pt,vh [2ll . |22j. The techniques 
for resumming both types of logarithms to all orders are well known and the fixed order 
perturbative and resummed calculations can be consistently matched at intermediate values 
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of the kinematic variables. 

We consider the process pp — > VH + X and present results from both the threshold 
resummation and the transverse momentum resummation of large logarithms separately for 
LHC energies. Since the final state particles are color-singlets, both types of resummation 



can be straight-forwardly adopted from results in the literature for t 



re Drell-Yan process [201 - 



24|. (We do not discuss the joint resummation of the logarithms [25j). Section HT1 contains 
a brief review of the resummation formalisms we apply. Details are relegated to several 
appendices. Section ITTT1 presents results for the total cross section including the resummation 
of threshold logarithms and a discussion of the theoretical uncertainties, while Sections IIV Al 
and IIV Bl contain some kinematic distributions resulting from the resummation of Pt.vh and 
threshold logarithms, respectively. Finally, Section |V] discusses the relevance of our results 
to searches at the LHC. 

II. RESUMMATION FORMALISM 

In this section we briefly review the transverse momentum and threshold resummation 
formalism that we utilize in deriving our numerical results. 

A. Transverse-Momentum Resummation 

The discussion of the transverse momentum resummation follows that of Grazzini et 
al. [^lj]. The hard scattering process under consideration is Higgs boson production in 
association with a vector boson in hadronic collisions 

AB V* + X -> VH + X, (1) 

where V = W ± , Z and X is the hadronic remnant of a collision. We apply the well known 
impact-parameter space (6-space) resummation (3, Q to the partonic cross section, 

J^. J~resum infinite 

a&VH _ da VH |_ a<y vH 

dMl H dp\ VH dMl H dp 2 T VH dM* H dp 2 TyH ' 

where Ptvh is the transverse momentum of the VH system and , Af a v ? $ contains the 

VH Pt,VH 

( M 2 \ 

resummation of the log a VH enhanced terms. Since all the logarithmically enhanced 

\ Pt,vh I 
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terms are factored into the resummed piece, the remaining contribution 
as Pt,vh —> and can be computed at fixed order in a s [2l| : 

davH 



(1(7 



~ finite 



VH 



dM VH d PT,VH 



is finite 



j a finite 

da J VH 



dM VH dp 2 TVH 



f-o 



dM VH dp 2 



■if-o 



J ~resum 
U0 VH 



dM 2 H dp 2 TVH 



(3) 



l VH"'fT,VH 

where the subscript f.o. refers to a fixed order expansion. In the low transverse momentum 
region, pt,vh "C Mvh, the resummed distribution is dominant, while in the high transverse 
momentum region, Pt,vh ~ Mvh, the perturbative expansion of the cross section dominates. 
Using Eq. flS]), the two regions can be consistently matched in the intermediate Pt,vh region, 
maintaining theoretical accuracy. 

To correctly account for momentum conservation, transverse momentum resummation is 
performed in impact-parameter space: 



M 2 VH - 



J " resum 
a<J VH 



M 2 VH 



db-J (bp TyH )W (b, M V H, s, fir, (if) 



(4) 



dM VH dp 2 TVH s j q 

where Jq{x) is the th order Bessel function and [i r , [if are the renormalization/factorization 
scales. By performing a Mellin transformation 1 it is possible to factor the terms that are 
finite and logarithmically enhanced as Pt,vh — >■ 0: 

M VH M VH M VH 



W% H {b,M V H,Vr,Vf) = H V N H (m vh 



x exp<^ G N ( a s (n r ),L 



fx r ' /if Q 
M YH M V h 



Q 



(5) 



where L = lnyQ 2 b 2 jb\ j with bo = 2exp(— 7^), H contains the finite hard scattering coeffi- 
cients, and G contains the process independent logarithmically enhanced terms. Hence, all 
the terms that are divergent as Pt,vh are exponentiated into the function GV, achiev- 
ing the all-orders resummation. The split between the finite and logarithmically enhanced 
terms is somewhat arbitrary; that is, a finite shift in the invariant mass Mvh can alter the 
separation: 

' Q 2 



log 



Pt,vh 



log 



Pt,vh 



+ log 



Q 2 



(6) 



The scale Q, termed the resummation scale, is introduced to parameterize this arbitrari- 
ness and is the same as that in Eq. (jSJ). To keep the separation between the finite and 



1 The Mellin transformation of a function h(z) is defined as /ijv = /q dzz N 1 h(z). 
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logarithmically enhanced terms meaningful, the scale Q has to be chosen to be close to 
M VH . 

As mentioned in the previous paragraph, all of the logarithmically enhanced contributions 
are contained in Gn- The divergent pieces can be reorganized such that Gn is written as 
an expansion that is order- by-order smaller by a s [2l| : 

G N [a S} L } , — — = Lg N (a s L) + > — g N [a B L, , — — , (7) 

V fir Q J t^2 Kl{) V Q J 

where g 7 ^ = for a s L = and Lg l N contains the leading log (LL) terms a™L n+1 , gj^ 
contains the next-to-leading log (NLL) terms oQL n , etc. Since the large logarithms are 
associated with collinear and soft divergences from real radiation, the functions g l N are only 
dependent on the initial state partons and are independent of the specific hard process under 
consideration. Explicit expressions for the LL and NLL terms needed for pp — > VH + X are 
given in Appendix A. 

The resummed distribution is valid in the low Pt,vh ^ Mvh region, while the perturba- 
tive expansion is valid in the high Pt,vh ~ Mvh region. However, as Qb approaches zero 
the logarithm L grows uncontrollably. As a result, the resummed distribution makes an 
unacceptably lar ge c ontribution to the high pt,vh region. This problem can be solved via 
the replacement 28] L — > L = \og^Q 2 b 2 /bl + l \ , such that L w L for Qb 3> 1 and L w 

for Qb <C 1. Hence, using L, the resummed contribution maintains the correct dependence 
on the large logarithms at low Pt,vh and does not make unwarranted contributions to the 
high Pt,vh region. This replacement has the added benefit of reproducing the correct fixed 
order cross section once the transverse momentum is integrated 2l|. 

The process-dependent function H is finite as Pt,vh 0. Hence, its Mellin transform 
H N does not contain any dependence on b and can be computed as an expansion in a s , 



tjVH 

/,„ M VH M VH M VH 



VH(n) ( M VH My H M VH 
fJ>r ' A*/ Q 



a (a s ,M VH ){l + J2(^) V 

^ 71=1 



(8) 



where do is the Born-level partonic cross section for qq' — > VH. At NLL accuracy, only the 
first hard coefficient H^ H ^ is needed. The value of this coefficient is given in Appendix A. 
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B. Threshold resummation 



In the original approach to threshold resummation lfj |l8| , the resummation is performed 



after taking the Mellin transformation of the hadronic cross section [29|, |30|. The Mellin- 
transformed hadronic cross section can then be factored into the product of the partonic 
cross section and the parton luminosity. The threshold logarithms for VH production are 
of the form ln(l — z), where z = My H /s, and are contained in the partonic cross section. 
After resummation, an inverse-Mellin transformation is performed to obtain the physical 
cross section. This leads to a new divergence due to the presence of the Landau pole in a s . 
Prescriptions for how to perform the inverse-Mellin transformation have been developed to 
remove this problem. The resummation of threshold logarithms for Drell-Yan production 
has been extensively studied 



31H34] 



More recently, techniques using soft-collinear effective theory (SCET) 35M38I] have been 
developed in which the resummation is performed in momentum space, obviating the need 

em of the Landau pole. In this paper, 



20| to the case of VH production. 



to go to Mellin space. This in turn removes the prob 
we will generalize the SCET resummation results of 

The leading singular terms at threshold in the hadronic differential cross section can be 
written as 

1 rlrr rl y /t- \ 

(9) 



rcr dM^ H 



— C(2, M VH ,fi f )C y-,Hf 



where r = My H /s with s the hadronic cm. energy-squared, £ is the parton luminosity, 



f 1 dx ( y \ , 

£(y,t*f) = J —fq( x ^f)fq'y-^fj +(g^q 



(10) 



and do is the Born level partonic cross section for qq' — > VH and is defined such that 
C(z, M VH , fif) = 5(1 — z) + 0(a s ). In the threshold region, z ~ 1, C(z, M VH , fif) can be 
factorized into a hard contribution and a soft contribution, 



C(z,M VH ,n f ) = n(M VH , v f )S{M VH {l - z),n f ). 



(11) 



The hard function %(Mvh, A*/) an d soft function S(M VH {1 — z),fj,f), evaluated at ///, are 
obtained by renormalization group running from the hard scale \lk ~ Mvh arid soft scale 
fi s ~ Mvh(1 ~ t), respectively, to sum the threshold logarithms to all orders in a s . 
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TABLE I: Approximation schemes for threshold resummation given a fixed order matched to a 
logarithmic approximation as in Eq. (|13p . 



Fixed order Log. Accuracy ~ a™L k T. 



cusp 



7 y ,7* C v , 



SDY 



LO NLL 2n - 1 < k < 2n 2-loop 1-loop tree-level 

NLO NNLL 2n - 3 < k < 2n 3-loop 2-loop 1-loop 
NNLO NNNLL 2n - 5 < k < 2n 4-loop 3-loop 2-loop 



The final result is found from that for the Drell-Yan process 20] 



c( z ,M VH ,fi f ) = \a 



VK 



-Mi 



VHi 



/j, h ) | U(M VH ,/j, h} ij, S} ij,f) — 



)l-2r/ 



xs DY In 



M VH {1 



fl 2 s Z 



-2y E ri 



T(2 V ) 



(12) 



where rj = 2ar(/x S) A 4 /), and Cy and sdy are the perturbatively calculable Wilson coefficient 
and soft Wilson loop coefficient, respectively. Eq. ( TTUj) with C given by Eq. (1T2|) is defined 
only for r/ > 0. For rj < 0, an analytic continuation is required. The analytic expressions for 
ar, Cy, §£)y and {7 which are necessary for our numerical calculations are given in Appendix 
B. 

Eq. (112j) is only valid in the threshold region z ~ 1. To obtain a formula valid for all 
values of z, we match the threshold-resummed result with the fixed-order result, 



da 



Here 



dM VH 



da 



matched 



da 



dM VH 



threshold resum 



da 



dM VH 



+ 



threshold f.o. 



da 



dM VH 



(13) 



f.o. 



is the result obtained using the threshold resummation formula of 



Eq. (H2X 



threshold resum, 
da 



VII 



is the fixed-order perturbative result and 



f.o. 



dM, 



da 
T~ 



is obtained 

threshold f.o. 

d singularity in C . The order 



from the fixed-order result by keeping only the leading thresho 
of the logarithmic approximation in the resummed result and the corresponding fixed-order 
results used in the matching of Eq. ( {TBI are summarized in Table [TJ 2 



2 The equivalence of the sub-leading logarithms between the SCET approach and the standard QCD Mellin 
transform approach has been studied in Ref . 39] . 
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III. SCALE DEPENDENCE OF THE CROSS SECTION 



In this section, we study the scale dependence of the total cross section for VH production 
at the LHC, beginning with the sensitivity of the resummed threshold distributions to the 
hard, soft, and factorization scales. Near the threshold, r = My H /s — > 1, the threshold 
logarithms are enhanced, leading to potentially large scale violations. The naive choice for 
the soft scale is /i s ~ My#(l — r). We follow the prescription of Ref. 20] to determine a 
sensible range of parameters for the soft scale. A low value of /z s is found empirically from 
the scale where the one-loop correction to Spy is minimal, 

M VH (1 - r) 



u (7) 



(14) 



2V1 + 100r 

Alternatively, an upper scale for the soft variation can be chosen as the value where the 
one-loop correction to sdy drops below 10%, 

M VH {1 - t) 



(ii) 



0.9 + 12r 



(15) 



.tun 



Empirically, the forms of ^ are insensitive to Mvh- Here and henceforth, we adopt the 
Higgs mass value 

Mh = 125 GeV. (16) 

We investigate the numerical effects of the scale variation by plotting the differential cross 
section of the threshold resummation of Eq. ( fl2|) and varying the soft, hard, and factorization 
scales. It is customary to measure the size of QCD corrections by a A-factor (A) typically 
defined as the ratio of a higher order cross section to the lowest order cross section: 

da „ da 



dMl H 



= A 



(17) 



LO 



where -tM, is a distribution defined at higher order in QCD. 

To study the scale variation arising from threshold resummation, we investigate the A- 
factor of Eq. ([17]) defined with 



da 



dM VH 



da 



dM VH 



(18) 



threshold— resum 

To isolate the effects of the scale variation due to threshold resummation from effects of the 
scale variation due to parton distribution functions (PDFs) and running a s , the A-factor is 



evaluated by using the NNLO MSTW20008 



401 ] PDF set and the 3-loop value of a s for all 



orders of the threshold resummed cross section and the LO cross section. Figure [T] shows 
the scale variation of this choice of i^-factor as a function of r at NLL between the dotted 
curves, NNLL between the dashed curves, and NNNLL between the solid curves for ZH 
production at Mzh — 1 TeV. The soft scale variation in pp — > ZH, with /i^ and /ij held 



constant, is shown in Fig. 1(a) The variation in the NLL result is significant, but the 
NNLL and NNNLL curves have little dependence on the soft scale, justifying the ad hoc 
choices of /zi 7 ' 77 ^. The K- factor grows rapidly as r increases, as expected. The sensitivity 



to the hard scale is shown in Fig. 1(b) with fixed /i s and [if. The hard scale is set by the 
invariant mass of the VH pair, and again we find that at NNLL and NNNLL, there is little 
dependence on /i^, showing excellent convergence of the perturbation series. Finally, we 



show the factorization scale dependence in Fig. 1(c) The factorization scale dependence is 
small even at NLL. 

We have also considered the scale dependence of the matched result for the total cross 
section. Analytic expressions for the LO and NLO fixed order results are found in Refs. 



llMl3|. and we use the computer code VH@NNLO for the fixed order NNLO results. The 



matched curves are found using the threshold resummation results of Eq. ( fl3|) . In Fig. 2, 
we use the MSTW2008 68% confidence level PDFs, and use LO PDFs for the LO and the 
NLL-LO matched curves, NLO PDFs for the NLO and NNLL-NLO matched curves, and 
NNLO PDFs for the NNLO and NNNLL-NNLO matched curves, and use 1, 2 and 3-loop 
evolution of a s respectively. We include the small contribution from the gg initial state in 
the ZH NNLO and NNNLL-NNLO matched curves. 



The results for ZH production at ^ = 8 TeV and = 14 TeV are shown in Figs. 2(a) 



and 2(b), respectively. We have chosen the central scale to be /i = M Z h- The top and 
bottom quark loops from the gg initial state contribute cr*'* loops = 0.06 pb at -y/s = 14 
TeV with fj,f = fi . This is the reason for the larger splitting between the NLO and NNLO 
curves than is seen in the WH results below. The fixed-order and matched curves have the 
renormalization/factorization scales set equal, fi r = fif. The matched and resummed curves 
have the hard scale, fit = 2Mvh, and the soft scale, fi s = §(/4 ) +^s H ^)- The NNNLL-NNLO 
matched curve is almost identical to the NNLO fixed order curve, and the resummation has 
little effect at this order. On the other hand, the NNLL-NLO matched curve increases the 
fixed order NLO result (at fif = /i ) by about 7%. 

The matched cross sections for WH production at y's = 8 TeV and 14 TeV are shown 
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H s dependence of pp— >ZH at M ZH = 1 TeV 



p., dependence of pp— >ZH at M =1 TeV 



(i) <ii) 



0.05 0.1 0.15 0.2 



T=M ZH 7s 



(a) 




0.3 0.35 



PL 
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M ZH <p h <2M ZH 
H f =M ZH ,p=(p s (I) +p s (n) )/2 



■■■■ NLL 
-- NNLL 
— NNNLL 



0.05 0.1 0.15 0.2 



T=M ZH 7s 



(b) 



[l dependence of pp— >ZH at M =1 TeV 



0.3 0.35 



0.5M ZH <^ f <2M ZH 
^= 2M ZH^ s =(H s <I,+ ^ s <n) ) /2 



NLL 

NNLL 

NNNLL 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 

T=M ZH 7s 

(c) 

FIG. 1: The (a) soft scale, (b) hard scale, and (c) factorization scale dependence of the threshold 
resummed cross section for pp —> ZH at NLL between the dotted lines dotted, NNLL between 
the dashed lines, and NNNLL between the solid lines normalized to the LO result (the if -factor is 
defined in Eqs. (|17p and (|18p ). The invariant mass Mzh is fixed at 1 TeV. 



in Figs. [2]^c) and[2^d). These figures show the sum of W + H and W~H production. As in 
the ZH case, the NNLO and NNNLL-NNLO matched results for WH production are quite 
close and show little scale variation. The NNLL resummation increases the NLO fixed order 
result by ~ 3%. 

The uncertainties in the ZH and WH cross sections from PDFs, renormalization and 
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pp->ZH+X, Vs=8 TeV, M H =125 GeV, Jl„=M z 



pp -> ZH+X, Vs=14 TeV, M H =125 GeV, H =M Z 




(a) 

pp^WH+X , Vs=8 TeV, M H =125 GeV, H„=M W 



(b) 

pp^WH+X, Vs=14 TeV, M H =125 GeV, (l„=M v 



Fixed-order - 
Matched 

0S L NNLL-NLO, matched NNNLL-NNLO, matched - 
^ I I 



NLO 



-4— 

NNLO 



NLL-LO, matched 



LO 



1.6 

a. l.s 

1 1.4 



NNLL-NLO, matched 
1 



NNNLL-NNLO, matched i . — 5— 

-- Fixed-order 

— Matched 




(c) (d) 

FIG. 2: Scale dependence of the fixed order (dashed) and threshold resummed matched (solid) 
cross sections for (a,b) ZH and (c,d) WH production at (a,c)-^/s = 8 TeV and (b,d) ^fs = 14 TeV. 
The NNLO and NNNLL-NNLO matched ZH results include the contribution from the gg initial 
state. 



factorization scale dependence, and the determination of a s have been investigated by the 
LHC Higgs Cross Section Working Group for the NNLO total cross section [7]. They find 
a total uncertainty at ^Js = 8 TeV of (9(4%) for WH and 0(5%) for ZH production for 
a 125 GeV Higgs boson. Our results show that including the resummation of threshold 
logarithms to NNNLL accuracy does not induce a ny further uncertainties. We note that 
Ref. {7] also includes the NLO electroweak effects 41], assuming complete factorization of 
the QCD and electroweak corrections. In the renormalization scheme, these corrections 
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reduce the total Higgs and vector boson associated rates by about (9(5%). 



IV. KINEMATIC DISTRIBUTIONS 



A. Transverse-Momentum Distributions 



We now give numerical results for the resummed transverse-momentum distributions. 
The distributions are computed at NLL-NLO accuracy with NLO MSTW2008 68% confi- 
dence level PDFs 



40| and the 2— loop evolution of a s using the formulae o 



Appendix A. 



The numerical results were found by modifying the program HqT2.0 [21| , |42| . |43|. The fac- 
torization and renormalization scales are set to the central values of fif = fi r = My + Mh- 
Also, the resummation scale is set equal to the invariant mass of the vector boson and Higgs 
pair, i.e., Q = M VH . 

Figures [3](a) andEJjb) show the transverse-momentum distribution for ZH and WH pro- 
duction, respectively, at = 14 TeV. The matched transverse-momentum distribution 
defined by Eqs. ([2]) and ([3]) (solid), resummed (dot-dash), fixed order expansion of the re- 
summed (dashed), and fixed-order perturbative (dotted) distributions are shown separately. 
As expected, the fixed-order expansion of the resummed and perturbative distributions are 
in good agreement. Hence, the finite piece, defined to be the difference between the pertur- 
bative distribution and fixed-order expansion of the resummed distribution as in Eq. (E]), is 
negligible at low transverse momentum and the matched distribution is dominated by the 
resummed contribution. The transverse-momentum distribution is peaked around 5 GeV 
for both WH and ZH production. 

For comparison, in Figs. [3]^c) and [3](d) we present the normalized matched transverse- 
momentum distributions for ZH and WH production, respectively, at both a/s = 8 TeV 
(dashed) and \fs = 14 TeV (solid). The position of the peak of the transverse distribution is 
not significantly different between the two LHC energies. However, the distribution at = 
14 TeV has a longer tail than at -y/i = 8 TeV. This can be understood by noting that higher 
transverse-momentum events correspond to higher partonic center-of-mass energies. Since 
events with higher partonic center-of-mass energies are more easily accessible at ^/i = 14 
TeV than at y/s = 8 TeV, we would expect there to be larger fraction of high Pt,vh events 
at y/s = 14 TeV than at y/s = 8 TeV. Hence, the transverse-momentum distribution has a 
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Q = M 7 



Vs = 14 TeV 
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Q=M WH M H =125GeV 

NLL-NLO 
MSTW2008 NLO 68%. 
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-■ Resummed 
-- Resummed f.o. 
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40 60 
P TWH [GeV] 

'(b) 



pp->W H + WH 



M-f =^ r = M W +M H 



M u = 125 GeV 

ri 




NLL-NLO 
MSTW2008 NLO 68% : 

-- Vs = 8TeV : 
— Vs = 14 TeV : 



40 60 
Ptwh [GeV] 

'(d) 



100 



FIG. 3: Transverse-momentum distributions for (a,c) ZH, and (b,d) WH production at the LHC. 
In (a) and (b), the matched distribution is shown with a solid line, the resummed distribution with 
a dot-dash line, the fixed-order expansion of the resummed distribution with a dashed line and 
the fixed-order perturbative distribution with a dotted line at y/s = 14 TeV. In (c) and (d), the 
normalized matched transverse-momentum distributions are shown for both y/s = 8 TeV (dashed) 
and y/s = 14 TeV (solid) LHC. 



longer tail for yfs = 14 TeV. 

Finally, we comment how the transverse-momentum resummation can effect the analysis 
of kinematical cuts on the signal cross section, particularly in relation to jet vetoes. At 
hadron machines, the VH production with Higgs decaying to bb has large QCD backgrounds. 
To reduce the backgrounds and effectively trigger on the signal, one usually considers leptonic 
decays of the vector boson. However, if the vector boson decay contains missing energy, 
W — > Iv or Z — > uV, semileptonic decays of tt can be a significant background. Since the tt 
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FIG. 4: Integrated matched transverse-momentum distributions normalized to the total cross 
section for both (a) ZH and (b) WH production. Results for y/s = 8 TeV and y/s = 14 TeV 
are shown with dashed and solid lines, respectively. 

background typically has more hard jets than the VH signal, a jet veto may be applied to 
suppress this background. We note that vetoing jets with a minimum transverse momentum 
can be approximated by placing an upper limit on the VH transverse momentum and, as 
can be seen in Figs. E^a) and[3](b), the perturbative calculation is unreliable in this regime. 
Hence, to fully account for the effects of a jet veto the soft-gluon resummation is needed. 
There has been much recent work on the systematic resummation of the large logarithms 



associated with jet vetoes 44j-|48 



To approximate the effect on the total cross section of a veto on jets with transverse 
momentum larger than p? vh we define 



T,VH) 



PT ' VH da 

dqr,VH-j , 

o dqT,VH 



(19) 



where da/dqxyH is the matched transverse-momentum distribution at NLL-NLO in Eqs. ([2]) 
and (El). Figure H] shows this cross section normalized to the total Pt,vh resummed and 
matched cross section as a function of Pt,vh for (a) ZH and (b) WH production for both 
y/s = 8 TeV (dashed) and y/s = 14 TeV (solid). As noted before in the discussion of 
Figs. [3](c) and |3]^d), at \fs = 14 TeV, there is expected to be a larger fraction of high 
transverse-momentum jets than at y/s = 8 TeV. Hence, <7(pt,vh)/& grows more slowly at 
y/s = 14 TeV than at y/s = 8 TeV. From the figures we see that the effects of a 20 (30) 
GeV Pt,vh cut decreases the NLO cross section by ~ 45% (~ 33%) and ~ 50% (~ 37%) at 
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y/s = 8 TeV and y/s = 14 TeV, respectively. 



B. Invariant-Mass Distributions 

In this section, we give numerical results for the invariant-mass distributions including 
threshold resummation and matching, using the analytic formulae of Appendix B. Since 
the distributions vary over many orders of magnitude, it is easier to see the effects in the 
.fT-factor, as defined in Eq. (1171) . Figures E(a) and^b) show the if-factor versus r at NNLL- 
NLO with a/s = 14 TeV for pp — > ZH+X and pp —> WH+X, respectively. The K-factor for 
the matched result of Eq. ffl3|) is shown with solid lines, the threshold-resummed contribution 
with dot-dashed lines, the fixed-order perturbative contribution with dashed lines, and the 
contribution from the leading threshold singularity of the fixed-order perturbative piece 
with dotted lines. Here we use MSTW2008 68% confidence level PDFs [40]. The scales are 
chosen to be \i$ = Mzh, Hh = 2Mzh and fi s = + /4 7 ) as in Section IITTI For the NLO 
fixed-order result, the leading threshold singularity of the NLO fixed-order result and the 
threshold-resummed result at NNLL, the NLO PDFs and 2-loop a s are used, whereas for 
the LO fixed-order denominator of the .ff -factor, we use the LO PDFs and 1-loop a s . As 
expected, the leading singularity and fixed-order results (the two lower curves) are close to 
each other, since the leading singularity dominates in the fixed-order result. On the other 
hand, the resummation effect is significant at high r, as seen by the large enhancement of 
the NNLL (the two upper curves) from the NLO result of ~ 20% for both ZH and WH at 
r = 0.3. 

The decrease of the X-factor at higher r values is due to the PDF effect. To see this, 
we artificially adopt the NLO MSTW2008 68% confidence level PDFs and 2-loop a s for the 
NLO fixed-order result, the leading threshold singularity of the NLO-fixed-order result and 
the threshold-resummed result at NNLL, as well as the LO denominator, and show the K- 
factors of these results with ^Js = 14 TeV for pp — > ZH + X and pp — > WH + X in Fig. |6f a) 
and Mjo) respectively. This is to isolate the effects of PDFs from a dynamical origin. The 
choice of scales is the same as in Fig. |5j We note that the monotonic increase of the K- 
factor distributions in Fig. |6] is drastically different from that in Fig. [5j This demonstrates 
the importance of a consistent choice of PDFs as in Fig. 

To examine the convergence of the perturbative series, we plot the i^-factors for the 
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FIG. 5: ET-factor distributions at y/s = 14 TeV for (a) ZH and (b) WH production. The NNLL- 
NLO matched result is shown with solid lines, the NNLL threshold resummed result with dot- 
dashed lines, the leading threshold singularity of the NLO fixed-order result with dotted lines, and 
the NLO fixed-order result with dashed lines. 

resummed results at NLL, NNLL and NNNLL with y/s = U TeV for pp -> ZH + X 
in Fig. using NNLO MSTW2008 68% confidence level PDFs and 3-loop a s for all the 
resummed results as well as the LO denominator. We see from Fig. [7] that the difference 
between NNLL and NNNLL is tiny (< 1%), confirming the excellent convergence of the 
perturbative series at this order especially after leaving out the PDF effect. 



V. CONCLUSIONS 

Given the exciting discovery of a Higgs-like particle at the LHC [l, 2], it becomes imper- 
ative to determine its properties. Thus its production rate at the LHC must be calculated 
as accurately as possible. Since the gauge boson-Higgs associated production is one of the 
channels that unambiguously probes the VVH coupling with V = or Z, it is of partic- 
ular interest. We combined the long- known fixed-order perturbative QCD calculations for 
VH production Jq] with soft-gluon resummation of both threshold logarithms and logarithms 
which are important at low transverse momentum of the VH pair. 

After a brief overview of the resummation formalism, we carried out detailed numerical 
analyses at the LHC for \fs = 8 TeV and 14 TeV. The overall corrections from NNLO fixed 
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FIG. 6: If -factor distributions at y/s = 14 TeV for (a) ZH and (b) WH production. The NNLL- 
NLO matched result is shown with solid lines, the NNLL threshold resummed result with dot- 
dashed lines, the leading threshold singularity of the NLO fixed-order result with dashed lines, and 
the NLO fixed-order result with dotted lines. The NLO PDFs and 2-loop a s are adopted for all 
the results as well as the LO denominator. 

pp^ZH+X, Vs=14 TeV, M H =125 GeV 




FIG. 7: pp — > ZH + X K factor distribution at yfs = 14 TeV for the threshold resummed piece at 
various orders of the logarithmic approximation, using the same PDFs for all curves. 

order calculations are sizable, increaing the LO rate by a factor as large as about 30% 
After implementing threshold resummation, the dependence of the cross section and various 
kinematic distributions on the soft and hard scales, as well as on the factorization scale is 
very weak, indicating the reliability of the calculations. The NNLL threshold resummed 
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total cross section increases the fixed-order NLO result by about 7%, while the NNNLL 
resummed result has little impact on the NNLO fixed order rate, demonstrating the excellent 
convergence of the perturbation series. 

The transverse-momentum spectrum of the VH system is calculated via soft and collinear 
gluon resummation. The distribution is peaked near 5 GeV and the spectrum is slightly 
harder at the center-of-mass energy of 14 TeV than at 8 TeV. Using the matched transverse- 
momentum distribution, we have also calculated the effect on the NLO cross section of 
placing an upper bound on the pr of the VH system. Since such an upper bound on the 
transverse momentum of the VH system limits the amount of transverse momentum a jet 
may carry in VH + X events, we expect the upper bound on the Pt,vh of the VH system 
to approximate a jet veto. 

As a final remark, our calculations can be easily extended to other electroweak pair pro- 
duction processes with the same color structures which arise via qq' annihilation at leading 
order, such as the EW gauge boson pairs and the Higgs pair production H°A°, H Q H^. A°H ± 

n 

and H + H- |49j. 



Appendix A: ptv Resummation 



In this appendix, we list the functions needed for the Pt vh resummation of Section III Al 



50 



511 ] . All formulae in this appendix can be found in Ref. 21[, but we include them for the 



convenience of the reader. First, the coefficients of the QCD beta function are normalized 
according to the expansion 

n=0 ^ 



(iln/x 2 



47T 



(20) 



At LL only the function g l N is needed and the Born level contribution arises only from qq' 



scattering 



2l| 



±A\\ A + ln(l - A) 



A = 



4tt 



a s (fi r )L, (3 



33 - 2n 



Al 



C f 



(21) 



andL = ln( 
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At NLL the functions g% and i?^ are needed 21]. 
M VH M VH , 



9N( a sL, 



ft, 



Q 



0i 



A, 



B 



q,N 



-^ln(l-A) i 

AA 1 



A 



+ ln(l - A) 



80 



2 / 153-19ry 

C F T67 7T 2 
IT 



1- A 



1^ + 27^ + ^ In ( 



M 2 VH 



Q 2 



(22) 



and rif is the number of light flavors. The anomalous dimensions are the Mellin 

transforms of the DGLAP splitting functions, P a b jo^l ]: 



lab,N 



n=l 



7T 



7aV = / dzz»^P ah {z) 
Jo 



(23) 



The process dependence arises through H 



VH(1) 
N 



53 



54| 



where, 



H 



VH(l) 
N,qq-^qg 



%Q,N lQ g -2 + 



rrVH(l) 
11 N,qq^qq 



li) 2(A + l)(A + 2) 

1 7T 2 \ 1 



(24) 



+ 



A(A + 1) 6 



+ -A 
2 



' 2 



.4 



VH 



+ 



2tt 5 



(26) 



Appendix B: Threshold Resummation 



In this appendix, we 



ist the functions needed for the threshold resummation of Section 



IHBl taken from Ref. 20j. All formulae in this appendix can be found in Ref. [20], but we 



include them for the convenience of the reader. 
The running kernel U is defined as 



U(M, n h ,ii s ,ii f ) 



exp [4S(fj. h ,fi s ) -2a 1 v(fJ, h ,fJ, s ) + 4a 7 ^>(/x s , ///)] , (27) 
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where a 7 is the anomalous exponent of 7 defined by 



(28) 



and jS is the Sudakov exponent 

5(z/, /i) 



c* s (aO 

a s (V) 



, -Tcusp(a) 



a s (i/) 



da' 

W) 



(29) 



The renormalization group equations, Eqs. fl28|) and (1291 . can be solved perturbatively. 
The anomalous dimensions are expanded as 



l{ot s 



7o^ + 7i 



47T 



+ 72 



4tJ 



+ 



The solutions to Eqs. (|28|) and (1291) are then 
a 7 (z/,/i) = 

+ 
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where r = a s {fi) / a s {p) 
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The cusp anomalous dimension is known to three-loops 55|, |56) . The coefficients are 



r = 4CV, 
ri = ac f 

To = 4CV 



r 2 



67 _ 7T 2 
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(33) 



The four-loop coefficient r 3 has not yet been calculated, so we use the Pade approximate 
T 3 = F^/Ti. The anomalous dimension 7^ can be obtained from the partial three-loop 
on-shell quark form factor [57]. The coefficients are 



7o y = -6CV, 
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(34) 



The final anomalous dimension, 7^, is known from the NNLO calculation of the Altarelli- 
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Parisi splitting function 



561 ] . The coefficients are 
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The other functions needed are the Wilson coefficient Cy and the soft function sdy- The 
Wilson coefficient Cy has the expansion, 

C v (-M 2 -ie,fi) = l + ^y (-L 2 + 3L-8 + ^Pj+C F (^y (C F H F + C A H A +T F n f H f ), 

(36) 

where L = ln(M 2 //n 2 ) — in, and 
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This agrees with the corresponding expression in [58l |. 
The soft function to two-loops is 

s DY (A n) = l + ^ (V + + (^) 2 (C F W F + C A W A + T F n f W f ), (38) 



22 



where 



2vr 2 it 4 
W F = 2£ 4 + ^ 2 + ^- 
o 1c 



w A 

Wf 



22 „ /134 2tt 2 . ., 
¥ (~9 ^ + 



9 9 27 



3 J V 27 
224 . 656 IOtt 2 8 



08 \ n 2428 67tt 2 
— + 28C 3 £ + + 

S3 I gl T 



+ -Cs- 
27 9 s 



7T 



22 



-Ca, 



(39) 



This again agrees with the moment-space expression in 581 ] 
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